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Abstract 

We present the exact and explicit solution of the principal chiral field in two dimensions 
for an infinitely large rank group manifold. The energy of the ground state is explicitly 
found for the external Noether's fields of an arbitrary magnitude. At small field we found an 
inverse logarithmic singularity in the ground state energy at the mass gap which indicates 
that at N = 00 the spectrum of the theory contains extended objects rather than pointlike 
particles. 
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1 Introduction 



The progress in the modern quantum field theory, especially in the solution of realistic prob- 
lems like QCD, quantum gravity and strings in physical dimensions, seems to be considerably 
slown down because of a very limited list of exactly solvable models of many interacting de- 
grees of freedom. Currently, this list includes the following relatively well explored systems: 

1. Quantum field theories and statistical mechanical spin models for the dimensions 
D < 2; 

2. Strings or 2d gravity coupled to matter with central charge c m < 1; 

3. Quantum gravity for d < 2. 

All the attempts to make some substantial progress on the way to realistic dimensions 
are usually confronted to the obstacles which every time seem to have a very similar math- 
ematical origin. It is clear that some principally new mathematical tools will be needed to 
attack in the future all these problems. In the absens of them we may try to use the old 
tools to get an insight to the physics of realistic dimensions for the string theory. 

One conceivable loophole would be to take some integrable two dimensional field theory 
of NxN matrix valued field. Its planar diagrams might describe fluctuating world sheets of 
a string propagating in the 1+1 dimensional space. 

Essentially, there exists only one candidate - the principal chiral field (PCF) on, say, 
SU (N) group manifold. It can be defined by the following action: 



where g is an NxN unitary matrix. Its large N solution has been anticipated for a long time 
to follow from its finite N solution || || Q . It was finally explicitly established in our paper 
[H]. We will follow the guidelines of this solution. 
The partition function of the model reads as: 



where the integration over g goes with the Haar measure for SU(N) at every point of the 
2D space. To test physics in this asymptotically free theory we have to introduce somehow 
an energy scale. Technically the simplest way to do it is to take the theory on a cylinder 
with the space compactified on the space interval [0, L] and to introduce the following twist 
of the field g: 




(1) 




(2) 



g(x,t) 



e HLX g{x,t)e H * 



(3) 
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where Hr(l) = diag(h\, h 2 — hi, fiN-i — — ^jv-i) is an element of the Cartan subal- 

gebra. 

It amounts to introducing in eq.(fl|) the covariant derivative 

D„g = d,g - i5, {H L g + gH R )/2 (4) 

instead of the usual derivative d^g. In what follows, we shall consider only the case H L = 
H R = H. 

In the hamiltonian language one adds to the hamiltonian of the PCF a term tr^HiQ^ + 
HrQr)/2 to the hamiltonian of the theory corresponding to the lagrangian (1). Here Ql = 
J cPxgdog' 1 and Qr = J d 2 xg~ x d§g are the Noether's left and right charges. 

The model is completely integrable for any N by means of the Bethe ansatz approach 
H or the bootstrap procedure for the physical S-matrix. The results show the presence of 
N — 1 physical particles in the spectrum of the theory obeying the masses: 

sm(f) 

where I — 1, N — 1 is the rank of a fundamental representation and m = mi is the mass 
of the vector particle. In the two-loop approximation it is 

1 — Att 

m = A — exp(— — ) (6) 
v Ao A 

where A is a cutoff. These results are in complete agreement with the renormalization group 
predictions for the asymptotically free theory. 

All particles are bound states of the vector particles. They form the multiplets of the fun- 
damental representations of SU(N) algebra (antisymmetric tensors corresponding to columns 
of the Young tableau). It follows from the last statement that hk is precisely the chemical 
potential of the k-th type of the physical particles. 

To demonstrate this last statement let us turn for a moment to the one- dimensional 
PCF and calculate iths partition function in the periodic time (0,T) with the same twisted 
boundary conditions. 

The functional integral representation of the partition function is 

Z(T, H L , H R ) = J [Dg(t)] H e mr So dt9 *9- l9 *9, (7) 

where [Dg(t)]R- is the Haar measure on the unitary group U(N), and the holonomies at the 
end points of the time interval are g(0) = c iHl and g(T) = e lHR (with e %Uh and e lHR the 
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eigenvalues of g(0) and g(T)). Due to the invariance of the Haar measure at each moment 
of time, we have, after the introduction of a new variable (connection) A(t) = ig~ l d t g(t) the 
following representation of Z(T, H L , H R ): 



Z(T,H ,H T ) = J DA(t)e- Ntr fo dtA W S^Te' ^ dtA(f) )e iH \ e iH 
Using the character expansion of the group 5-function: 

5(U,U') = Y J XR{U)x R {Uy (9) 

R 

and the fact that the character is the trace of the matrix element in a given representation 
R: xr(T exp(2 / AM)) = tr R [T exp(i / A a T%dt)\ (rg is the ath generator of U(N) in the Rth 
irreducible representation), we arrive after a simple gaussian integration in A(t) (independent 
at each moment t) at the conclusion that 

Z(T,H L ,H R ) =Y,z- TC2iR) XR{e^)xR{e lHR ) (10) 

R 

where C^-R) is the second Casimir of the group U(N) in the representation R. This is 
also equivalent to the partition function of two dimensional QCD on the cylinder with the 
boundary holonomies specified at either end of the cylinder by the two distributions e iHh 
and e lHR . 

Let us now recall the Weil formula for the character: 

iH _ det fcJ (e^~ fc +^) 
)- , (11) 

where A(x) = l\i<j( x i ~ x j) is the Van-der-Monde determinant. 



Let us consider one of two characters entering the expression (1C). Due to the antisym- 
metry of the characters in n^'s we can retain in one of them only a diagonal term from the 
Weil determinant under the sum over irreps (in other words, over integers n^): 

exp(i H kUk) = exp(i ^2(h k - h k ^ x )n k ) = exp(-z h k(n k +i - n k )) (12) 

k k k 

Since n k +i — n k is the number of columns of length k in the Young tableau of the irrep R 
(corresponding to the antisymmetric group of k indices from the tensor describing the whole 
irrep R), we conclude that h k is the chemical potential of k- vectors - the k-th fundamental 
representations which constitute the whole irrep R. In the 2D PCF these k-vectors will 
correspond to the physical particles of the type k. 
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2 Ground state energy and beta-function of the PCF 
at large N 



Let us turn back to two the two-dimensional PCF, formulate our main result and then discuss 
it in details. 

Let us make a special choice of h^s which will be technically the simples one and quite 
suitable for the limit N = oo: 

= (13) 

Since every k-th type of particles is excited only if hk exceeds rrik, we will have for this 
choice of hk's no particles excited when h < m, and all particles excited on equal footing 
when h > in. 

We show that the energy of the ground state is expressed in terms of modified Bessel 
functions: 

f(h) = -L(8(h) - 8(0)) = —^B 2 Ix(B)Ki(B) (14) 



N 2 

where the parameter B is defined through 

171 

- = BK\{B) (15) 

The distribution of rapidities of physical particles will obey the simple semi-circle law with 
the support B. The parameter B defines the value of rapidity corresponding to the Fermi 
momentum of the fused particles. We shall see, that B gives the most natural definition of 
the renormalized (running) coupling constant : 

X(h) = ^ (16) 

With the definition flUf ) of the running charge one can find from eq. (|T5| ) the exact 
beta-function: 

, d - 4tt dB A K X {B) 



or 



^-t^r-l~ x tM^ n (18) 



where 



&o = l, 6i = 4, 6 2 = -8,6 3 = 64, 6 4 = -5 2 -2 5 , 
& 5 = 13 ■ 2 10 , 6 6 = -1073-2 8 , ,6 7 = 103-2 16 ' ' 



b n ~ -(-l)" v /8/( 7 rn)(4n/e) n , n -> oo (20) 
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3 Physical consequences of the exact solutionrweek and 
strong coupling limits 



Let us first look for the week coupling regime and find the link with the results of the 
standard renormalized perturbation theory. In the asymptotically free theory week coupling 
means the presence of a big energy scale, namely, big h/m. From eq. ([15]) we conclude that 
it corresponds to B —>■ oo. It follows from the large B asymptotics of McDonald function 

^ (l-^ + 0(l/5 2 )) (21) 



m V7Tv / S v %B 

Solving for B we obtain: 



B = \ n — + -\n\n — + -\n- + 0(-^j-) (22) 
m 2 m 2 2 v lnA ; K J 

m 



Using the large B asymptotics h(B) we finally find from (HP and (|2 

lQnf(h) = -h 2 B + O& = -/i 2 (ln — + ilnln — + iln^ + 0( * )) (23) 

B v m 2 m 2 2 m — ' 

m 

This result reproduces correctly one- and two-loop terms of the perturbation theory as well 
as the universal non-perturbative constant | In |. We also see from the comparison of eq. (|22|) 



with the standart 2-loop perturbation theory that our definition (|16D of the running coupling 
was justified. 

More than that: we can use the known formulae for the 1/B expansion of the product 
Ii(B)Ki(B) to get the explicite coefficients of the renormalized week coupling expansion: 

/C)A 2 = -i(l-£G 2 „(^) 2 ") (24) 



where 



2n+l((2n)!) 3 /=r^2n 

L/Ot? = ~, — ,s , „ >r)^oo 2 \/im — 25 

2n-l n! 4 8 2 " /v V e ; V ; 



Note that the expansion goes only in even powers. 

As we see, in spite of the fact that every coefficient represents a sum over renormalized 
planar graphs, it grows factorially with the order. Most probably this happens because of 
the renormalons (some subsequence of logarithmically divergent graphs) giving the main 
factorial contribution in each order noticed long time ago by 'tHooft. This means that we 
have an exponential number of graphs in each order but some of them give (2n)! contribution 
after the momenta integration. More than that: the series is a non-signchanging one and 
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thus non-Borel summable. Nevertheless, the free energy perfectly exists for any finite A. 
These phenomena seem to be imminent for any asymptotically free field theory. It sheds 
some doubts on the possibility to interpret the standard planar Feynman diagrams in this 
theory as the fluctuating world sheets of some string. The point is that in all examples 
of the solvable theories of non-critical strings with the central charge of the matter c < 1 
the expansions in the cosmological constant (conjugated to the invariant area of the world 
sheet, or, in the discretized version of the string theory, to the number of vertices of the 
corresponding random (Feynman) graph) have usually a finite radius of convergency, which 
means that near the critical value of the cosmological coupling the graphs are big. In the 
present case we have a zero radius of convergency with respect to the renormalized coupling 
(because of renormalons) , which prevents the usual procedure of the thermodynamical limit 
of a big world sheet. 

The hope to interpret the model in terms of a new string theory gets revived if one 
considers an opposite, nonperturbative limit of small energies (of the order of the mass gap): 
A = h/m — 1 — >• 0). This corresponds to B — > 0. 

At small B asymptotics of the Bessel functions give a singular behaviour on the threshold: 

f(h) ~ -(m/2ir) 2 j^- A^O (26) 

It differs drastically from the threshold behaviour for a finite N theory of massive particles, 
where we would have 3/2 law (see e.g. [[J): 

f N (h) ~ -m 2 (A) 3 / 2 (27) 

The reason for that is that we fixed the mass scale at mjv/2 at large N, and therefore the 
mass spectrum became continuous (we have N-l different masses on the finite interval). The 
simpler limit would correspond to fixing m\. In this case rrik = km\ and we would not have 
any bound states. 

It is interesting to compare the result (|26|) with the c = 1 matrix model where the 
ground state energy behaves in a similar way with respect to the cosmological constant A: 

EW -\£0h' (28) 

The mechanism by which the inverse logarithmic behaviour with respect to the cos- 
mological constant occurs also requires a parametrization through the fermi level of the 
corresponding fermions (whose coordinates are the eigenvalues of a hermitean matrix field). 
The Fermi level plays the role of a "hidden" parameter of the problem and the eigenvalues 
give rise to an extra (Liouville) degree of freedom of the theory. 
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4 Sketch of the solution 



Let us now sketch out the major steps of the solution of the PCF. 

The easiest way to start with is the construction of the physical S-matrix for the (vector x vector) 
particles [|, |]. Following the general guidelines of the bootstrap method one arrives to the 
result : 

S = X(6)S(6) ® S(6). Here 9 is the rapidity of a massive relativistic particle (p = 
mcosh^p 1 = msinh#) and X(6) is the CDD-ambiguity factor which cannot be determined 
by the factorization, unitary and crossing symmetry conditions. The SU(N) unitary, crossing 
invariant, factorized S-matrix of vector particles is well known. It is 

where is the projection operator onto symmetric (antisymmetric) states. 

Y(\ - -i-lrC-i- + -M 
u{6) = 2 r T a ( 3 °) 

Finally, the CDD factor X is chosen to cancel all double zeros and double poles on the 
physical sheet < ImO < n: 

yftft - sinh <f + I] r,n 

- ShtFf) (31) 

This is the 5-matrix of the vector particles. It has a pole on the physical sheet at 9^ = 2iri/N 
in the antisymmetric channel. It corresponds to the first bound state ( the second rank 
antisymmetric tensor) with a mass m 2 = msin(27r/iV)/ sin^/iV). The S'-matrix of these 
particles can be also found by tensoring the vector S'-matrix (the fusion procedure). It also 
has a pole in the antisymmetric channel, and so on. In this way the whole mass spectrum 
(H) can be generated. 

To find the thermodynamical properties of the ground state we have to use the Bethe 
ansatz procedure. The idea is to put many vector particles into a periodic box of the size 
L. In the thermodynamical limit they will self-organize into higher bound states, and to 
find the ground state energy we have to calculate the densities of physical particles Pk(0) for 
every type k — 1, 2, N — 1. 

The next step (which can be justified within the thorough consideration of the Bethe 
procedure) consists from throwing away the isotopic index structure of the S-matrix and 
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using its scalar part S(9) = u 2 (9)X(9) to write down the phase ballance eq. for a particle 
moving around the box: 

M 

exp(imL smb. 9 a ) = exp{i(p{9 a — 9 p)) (32) 

where exp{i(j)(9)) = u 2 (9)X(9) . 

To restore the whole SU(N) x SU (N) structure of the physical states of the model we have 
to consider the solutions of this eq. for complex values of #'s. It is known that in the large 
L limit the vector particles form the "strings" of complex rapidities: 9 r '^ — > + 2mi/N, 
where 9 l is a rapidity of the Z-th particle and r is an integer running between —1/2 and 
1/2. Substituting this into the Eq.(^) and multiplying equations over r we shall obtain the 
equations for the rapidities of the state which contains Mi particles of the kind /. Taking the 
logarithm of both sides the Eq. ( p2|) we obtain 

N-l Mn 

Lrm sinh0« = 27rJ«+X; £ toJfi® - 9f) (33) 

n=l a=l,^/3 

where <pi n (9) = J2\ r \<i/2,\r'\<n/2 + 2rm/N + 2r'm/N) is the scattering phase of the Z-th 
and the n-th particles, and integers J are the quantum numbers of the states. The energy 
of this state is obviously 

i N-l Ni 

E = -J2 m iY. coshe « ( 34 ) 

L 1=1 a=l 

One can note here that the ground state of the system can be described by a big Young 
tableau where every column of a length k corresponds to a physical particle of the type k 
which transforms as k-th fundamental representation of SU(N). 

The next step is to find rapidities to minimize the energy (|34]) in the thermodynamic 
limit Mi/L = ni, while L — ► oo. We assume that in the ground state #'s are distributed 
smoothly between —B\ and B\ with a distribution functions pi(9) Than eq.(^) implies the 
spectral equations 

^-mi cosh# = £ 12 R in{Q - 0')pn(9')d9' (35) 

where Ri n {9) = 5i n — ^d<pi n (9) / d9 . The Fermi rapidities Bi are determined by the number 
of particles in the / th representation: f^L pi(9)d9 = n/. The energy of the state is then 

E = £ f Bl rrn coshtf p t (9)d9 (36) 
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An explicit form of the scattering kernel Ri n was found in 0. Its Fourier transformation is 

sinh - ±j) sinh fe) 

R*n(o>) = 2 V . »' } (37) 

sinh TTUJ 

at I > n and R\ n = R n i. 

Now we pass to the large N solution (found in [IJ) of this rather complicated system of 
integral equations. 

At large N we can consider a particular distribution of fields hi which creates all different 
particles on equal footing, namely one which follows the spectrum of masses hi — 
(h/m)mi . This field creates Mi = Af(mi/m) particles in the l-th representation (the most 
representative Young tableau). In this case all Fermi momenta are equal: Bi = B and 
pi = jj{mi/m)p. Then the spectral equations fl35D can be easely diagonalized. They reflect 
the structure of the Cartan matrix and moreover have the same eigenvectors: 



N-l 

E 

Ln=l 



E X^RlnM^ = RP>W 



it^H = — V - } ul at\2 ( 38 ) 
V ; 7T J? uj 2 + (p + rN) 2 v ; 



where xf^ = sin^, p = 1,2,...,N — 1 (x^ is t ne Perron-Frobenius mode) is the 

orthogonal set of eigenfunctions of the Cartan matrix Ci n = 25i n — Si jTl+ i — 5j+i, n - 

Let us note that the expression (|38|) for the kernel of the integral equation reminds the 
propagator of the free motion of a particle on the space consisting from the physical space 
(the time plus the 1-dimensional physical space described by rapidities) plus the discrete 
periodized space represented by the Dynkin diagram. In this sense it looks like the effective 
space of our model becomes three dimensional, the phenomenon familiar to us from the 
experience with the c=l bosonic string (where it effectively passes from ID to 2D). 

Then the density p obeys the equation: 

- ( B R w (6 - 9')p(9')d9' = — cosh 6. (39) 
N J-B 2n 

Further simplifications occur in the large iV limit 

«<"M^J^ (40, 
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Now the density p may be found in a closed form. To see this, let us apply the operator 

a 2 
de 2 



+ 1) on both sides of the equation. As a result we obtain an integral equation with 



the Cauchy kernel (9 — 9') . This equation is solvable: 

Tfi 

p(6) = , (41) 

where K (B) is the Bessel function. Note that, in the large N limit R^(tu) vanishes at large 
u, whereas at finite iV it approaches 1 (see eq.(28)). This implies a singular behaviour of 
p{6) at the Fermi point ±B. As a result the physics on the threshold h ~ m will be changed 
drastically. 

The value of the Fermi rapidities as a function of number of particles can now be obtained 
from 

/B TTHl 
. B ^=*m (42) 

In turn the energy of the state with a given number of particles is 

And finally, the field h = —2^ 2 /N 2 dE/dn which corresponds to a given number of particles 
and the energy as a function of the field £ = E — J2i hini = E — N 2 /(2tc 2 ) hn are given by 
the formulae ( |15] ) and (PI). 

Let us note that in the large N limit any virtual and real processes involve all particles, 
since minimal energies are greater then a minimal separation between masses. A reasonable 
external field hi ~ mi excites all of them on equal footing and leads to collective effects. 



5 Conclusions 

Let us now list the main lessons which can be drawn out from the solution: 

1. Renormalized planar graphs do not model the world sheets of any string (they have 
the zero radius of convergency, due to renormalons) . 

2. A possible "stringy" behaviour with the 2 dimensional physical target space may occur 
only in the non-perturbative strong coupling regime (near the threshold {h ~ m). In this do- 
main the physics looks similar to the c = 1 bosonic string described by the matrix quantum 
mechanics 0. At least, the ground state energy contains the similar inverse logarithmic be- 
haviour with the scale paramenter, and the fermionic spectrum of excitations which becomes 
classical in the large N limit, is also common for two models. 
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3. The Bethe Ansatz solution of the PCF can be effectively described as having (1+1 + 1) 
dimensions: (time x space xDynkin diagram), according to the form of the integral equations 
for densities of particles. The picture is similar to (time <8> Liouville (eigenvalue) mode) 
forming a (1+1) dimensional effective space in the c=l bosonic string. 

Among the most obvious questions which are left we list the following ones: 

1. Elaboration of the 1/N expansion and of the double scaling limit around the threshold 
h ~ m, where a possibility for the "stringy" behaviour exists. 

2. One has to understand the role of the 3-rd dimension (Dynkin diagram) in the model. 
May be, the theory is "almost free" in the (1+1+1) dimensional effective space. 

3. Main challenge: to solve the model by simpler (matrix model) methods directly in 
terms of the original chiral field g(x). The character expansion used for the ID principal 
chiral field (equivalent to the heat kernel on the group SU(N) manifold) |7|, |J might be very 
useful for it. 
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